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Abstract 
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cone QCD sum rule approach. Comparison of our results on form factors with the 
existing results and lattice calculations are presented. 
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1 Introduction 



It is well known that the structure of the nucleon is parametrized in terms of form factors. 
The electromagnetic form factors of nucleon are measured in a wide range of momentum 
transfer squared q 2 (see for example [1] and the references therein). But, in contrast to 
the electromagnetic case, the form factors Ga(q 2 ), Gp(q 2 ) and Griq 2 ) due to isovector 
axial-vector current are not known. The nucleon matrix elements of axial-vector currents 
at q 2 = are determined by the axial-vector coupling constants such as gA (isovector), g\ 
(octet), g^ (isoscalar) and g\ (flavor singlet). Among these only nucleon isovector coupling 
constant gA is well known which is measured from the neutron /3-decay. Knowledge of any 
three of these coupling constants determines the quark spin content of the nucleon. For 
this reason study of these constants receives great interest. 

Using the Lorentz covariance, the matrix element of the isovector axial- vector current 
between initial and final states is parametrized as 

(N{p') \A,\N(p))=u(p') 

where = uj^jsu — dj^d, q = p — p', tun is the nucleon mass, and Ga, Gp and Gt are 
the axial, induced pseudoscalar and induced tensorial form factors, respectively, induced 
by the isovector and axial-vector currents. The matrix element of isoscalar axial-vector 
current A^ = vrf^^u + dj^^d between nucleon states is determined similar to (CQ) with 
the following replacements: Ga — > G% Gp — > Gp and Gt — > G T (here the superscript S 
means isoscalar). The form factors Gt and G T both vanish as a result of the exact isospin 
symmetry and G-parity invariance of the strong interaction. 

The aim of the present work is the calculation of the form factors Ga{G^) and Gp(Gp) in 
the framework of light cone QCD sum rules (LCSR) [3J H] approach, using the most general 
form of nucleon interpolating current. LCSR is based on the operator product expansion 
over twist of the operators near the light cone. This method combines the standard sum 
rules technique [5] with the parton distribution amplitudes describing the hard exclusive 
processes. This method is widely applied to the problems in the meson sector (see [1]). 
Lately, the electromagnetic form factors of nucleon [B] , the scalar form factor of nucleon [7] 
and the weak A& — > piv [S] are investigated in the baryonic sector in the framework of this 
method. Note that higher twist amplitudes of the nucleon are calculated in [9]. 

The outline of this work is as follows. In section 2 the basic ingredients for calculating 
the form factors due to axial-vector current are introduced and sum rules for the form 
factors are constructed. In section 3 we present our numerical results together with the 
concluding remarks. 

2 Sum rules for the form factors of the nucleon due 
to isovector and isoscalar currents 

In this section, we construct sum rules for the form factors of the nucleon due to the 
isovector and isoscalar axial-vector currents. For this purpose we start by considering the 
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polarization operator, which is the basic object of the sum rule approach, 



n>,g) = i / d 4 xeW0 T{ v (0)M s \x)} N(p) 



(2) 



where 



Ajf) = uj^u =F d'j^d 



with the upper (lower) sign corresponding to the isovector (isoscalar) axial- vector current, 
and 7] is an interpolating current with nucleon quantum numbers. The nucleon interpolating 
currents without the derivative terms, and with nucleon quantum numbers, can be written 
as [10] 



rji{ x ) 



2e abc Y, (u aT (x)CA[d b (x))A f 2 u c (x) 



r h {x) = -e afec [(M aT (x)C^(x)) 75 ^ c (x) - (M aT (x)C^(x)) 75 KW 



(3) 
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where A\ = I , A\ = 75 , A\ = 75 , A\ = /3; C is the charge conjugation operator; a, b, c 
are the color indices, and z is a light-like vector with z 2 = 0. The choice (5 = — 1 in Eq. 
([3]) corresponds to the Ioffe current Note that 772 current is modified in the following 
way [9] 
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The axial form factor and induced pseudoscalar form factor of the nucleon are calculated in 
[12] in the framework of LCSR using the current 773(2;). However the current 7/3 (x) couples 
to both spin 1/2 and 3/2 baryons. Therefore, for a reliable determination of the form 
factors the unwanted contributions coming from spin 3/2 states should be eliminated. But 
this elimination is not done in [12] . 

Note that the form factors Ga{Q 2 ) and Gp(Q 2 ) induced by the isovector axial current 
are calculated in [13] in LCSR using the Ioffe current. 

In the present work we calculate the form factors induced by isovector and isoscalar 
axial-vector currents using the general form of the interpolating 771 (x). 

We start by calculating the correlation function ([1]) from the QCD side. At large Eu- 
clidean momenta p' 2 = (p — q) 2 and q 2 = —Q 2 the correlation function can be calculated 
perturbatively. Using the expressions of interpolating currents we get the following result 
for the correlator: 
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where S q ( — x) is the light quark propagator and its light cone expanded expression is [13] 
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The terms proportional to G^ v in Eq. <£jty give contribution to four- and five-particle 
nucleon distribution functions, and these amplitudes are expected to be small [HI [15] 
which will be neglected in further analysis. Only the first term in Eq. (J7|) survives since 
Borel transformation kills the second term. 

It follows from Eq. that for the calculation of I1 M we need to know the matrix element 
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Ua(aix)up(a2x)d^(a 3 x) N(p)j. This matrix element of nonlocal operator is defined 
in terms of the nucleon distribution amplitudes (DAs), and their explicit expressions are 
presented in M, Mi EH HE] . 

Using the explicit expressions of the nucleon DAs and performing integration over x and 
selecting the structures 5^7^75 for Ga(G^) and q^^ for Gp(Gp), we get 
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2(1- (3)F 15 (t 3 ) + (I +(3) F 1G (t 3 ) 
where upper (lower) sign corresponds to isovector (isoscalar) axial- vector current, and 
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Those functions appearing in Eqs. (JED and (jHJ) are determined as 
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whose explicit expressions are given in [T3] . 

Physical part of the correlation function is obtained by inserting a complete set of states 
between the currents in Eq. (JTJ) with the same quantum numbers of the current rj(x). After 
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isolating the pole term of the nucleoli state, the correlator function ([I]) can be written as 
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where p' = p — q and q is the momentum carried by the axial-vector current. The second 
term in ([TO]) describes the higher states and continuum contribution, and h is complete set 
of the hadrons with the quantum numbers of the ground state nucleon. 

Contribution of higher states to the physical part of the sum rules are taken into ac- 
count using quark-hadron duality, i.e., spectral density for higher states is equal to the 
perturbative spectral density starting from s > s , where s is the continuum threshold. 

The matrix elements entering to Eq. (fTOl are defined as 
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Substituting Eqs. ( fTTj) and ( fT2|) into Eq. ( fTUj) and selecting the structures ^7^75 and q^js 
and performing Borel transformation with respect to the variable iq — p) 2 , which suppresses 
the continuum and higher state contributions, we get the above-mentioned form factors: 
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In performing Borel transformation, we use the following substitution rules (see for 
example [E3 CEBJ and [IS]) 



dx 
dx 



p{x) 



(q — xp) 2 
(q — xp) 4 



dx 



x 



p(x)e 



-s/M 2 



dx 

M 2 J 1c 2 



p(x)e 



-s/M 2 



+ 



Q 2 + x\m 2 N 



-so/M 2 



(15) 



where 



s(x) 



[1 - x)m N + 



1 — x 



x 



Q 2 



and xo is the solution of the quadratic equation for s = so, i.e., 

1 



x 



2m 2 N 



{Q 2 + s - m 2 N ) 2 + Am 2 N Q 2 - (Q 2 + s - m 



NJ 



and Q 5 



3 Results and discussion 

Now we are ready to examine the sum rules for the form factors. It follows from expressions 
of the sum rules for form factors that DAs are the main input parameters of LCSR. The 
complete list of all DAs which enter to the sum rules for form factors can be found in [IB] . 

These DAs contain eight hadronic parameters f^, Xi, A 2 , V^, A", ff, f$ and /" which 
should be determined in the framework of various models. In further numerical calculations 
we consider three different sets of these parameters: 
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• QCD sum rules based DAs, in which the parameters are determined from QCD sum 
rules (set 1) having the values A\ = 0.38 q= 0.15, Vf = 0.23 q= 0.03, ff = 0.40 q= 0.05, 
ff = 0.22 =f 0.05, fi = 0.07 q= 0.05. 

• A model for nucleon DAs (set 2) in which the above-mentioned parameters are chosen 
in such a way that the form of nucleon DAs describe well the existing experimental 
data on nucleon form factors (see p2]), whose numerical values are A\ = 1/4, Vf = 
13/42 and the values of ff, ff and /" are the same as in set 1. 

• Asymptotic forms of DAs of all twists (set 3) in which the values of the parameters 
are set to A" = 0, Vf = 1/3, ff = 3/10, ff = 4/15, ft = 1/10. 

Note that the values of f^, Ai and A2 for these three sets are the same, i.e., /jy = 
(5.0 =F 0.5) x 10^ 3 GeV 2 , X 1 = -(2.7 =F 0.9) x 10~ 2 GeV 2 and A 2 = (5.4 =f 1.9) x 10~ 2 GeV 2 . 

The values of the non-perturbative parameters entering to DAs at [i = 1 GeV scale are 
given in [9] and |13j . Sum rules for the form factors involve three unphysical parameters, 
namely, continuum threshold sq, Borel parameter M 2 and the parameter (3 in the interpo- 
lating current of nucleon. Of course, if we can perform OPE up to infinite order, the result 
must be independent of these parameters. But we truncated OPE in the finite order, and 
hence, there appears dependence of sum rules on these auxiliary parameters. However, any 
physical quantity can not depend on unphysical parameters. Therefore, in the first hand, 
we should find the appropriate regions of unphysical parameters where form factors are 
independent on them. From an analysis of mass sum rules it follows that, when continuum 
threshold sq lies in the region 2 GeV 2 <s < 2.5 GeV 2 , the prediction of sum rules on the 
mass of the baryons are practically independent of sq. It is this region of Sq which we will 
use in our numerical calculations. 

Having determined Sq, we next try to find a region (the so called working region) of M 2 
where the above-mentioned form factors are independent of M 2 at fixed values of so and 
(3. We study the dependence of the form factors induced by the axial-vector current on M 2 
at fixed values of Q 2 and (3 at so = 2.0 GeV 2 and so = 2.5 GeV 2 , for three different sets of 
the wave functions. We obtain that the results are almost the same for both choices of so- 

An upper bound for the Borel parameter M 2 is determined by requiring that the con- 
tribution of continuum be less compared to the continuum subtracted sum rules. Lower 
limit is determined from the condition that the contribution of the term with highest power 
of 1/M 2 is less, say 30%, compared to the higher powers of M 2 term contribution. Using 
these constraints we found that the working region of M 2 is 1 GeV 2 < M 2 < 2.5 GeV 2 , 
and we obtain that the results are rather stable with respect to the variations of M 2 , when 
M 2 lies in the above-mentioned working region. It is also observed that the results of for 
the form factors depend also on the parameter (3. But, as has already been mentioned, the 
parameter (3 is an auxiliary quantity and therefore the form factors must be independent 
of it. As a result of this argument, we need to find a region for (3, where the results for the 
form factors are independent of its value. Analyses of mass sum rules [19] and meson-octet 
baryon couplings [20] lead to the result (3 < —1.3 and (3 > 3.3. In our numerical calculations 
we will use these bounds for (3. It should be noted here that the Monte Carlo analysis of 
mass sum rules for baryons [21] predicts an optimal value (3 ~ —1.2, which is close to our 
choice of lower bound. 
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The correlation functions can be calculated in QCD for sufficiently large negative values 
of Q 2 and (p — q) 2 using OPE. The form factors can reliably be determined at the range 
Q 2 > 2 GeV 2 . Our approach is not applicable for smaller values of Q 2 . For this reason, 
the form factors are evaluated in the range Q 2 > 2 GeV 2 . We study the dependence of the 
form factors on Q 2 at fixed values of (3, lying within the above-mentioned working region, 
at fixed values of M 2 = 2 GeV 2 and s = 2.25 GeV 2 , for the three sets of DAs. 

From numerical analysis we obtain the following results: 

• The form factors Ga{Q 2 ) and Gp(Q 2 ) exhibit practically the same Q 2 dependence for 
all three sets of the nucleon DAs. 

• The values of the form factor Ga{Q 2 ) almost coincide at negative values of (3 for the 
second and third set of DAs, while it differs about 50% at (3 = —1.4 and 15% at 
(3 = —5, for the first set of DAs. For the positive values of (3, Ga(Q 2 ) practically 
coincide for all three sets of DAs. 

• The values of Gp(Q 2 ) are very close to each other for all three sets of DAs. 

• The values of G S A {Q 2 ) are very close to each other for the second and third sets of 
DAs, but it is larger about 50% for the first set of DAs at f3 — —1.4. 

• The difference between the values of Gp(Q 2 ) is very small for all three sets. 

• We see that all form factors are negative (positive) at positive (negative) values of (3. 

For illustration, in Figs. (l)-(4) we present the dependencies of the form factors Ga(Q 2 ), 
Gp(Q 2 ), G S A {Q 2 ) and Gp(Q 2 ) on Q 2 at several fixed values of (3, including the results for 
(3 = — 1 case. 

As has already been noted, the axial form factor Ga{Q 2 ) is considered in [13], however 
the analytical results that are presented in that work and in ours should be different, since 
the considered structures are different. However, we see that our numerical results on 
Ga(Q 2 ) are close to that obtained in [13j at j3 = — 1. 

Note that for the analysis of the existing data from neutrino scattering experiments, 
form factor Ga{Q 2 ) is usually parametrized in the dipole form 

Gf{Q 2 ) = °* j, 

(1 + Q 2 /m 2 A ) 2 

where gA = 1.2695 =F 0.0029, which is determined from {3 decay [2]. The global average 
for rriA extracted from neutrino scattering experiment is predicted to have the value rriA = 
(1.026 =f 0.021) GeV in (22], while the value m A = (1.20 q= 0.12) GeV announced by the 
K2K collaboration is slightly larger |23j . 

In Figs. (5)-(10) we present the LCSR prediction for the axial form factor Ga{Q 2 ) nor- 
malized to Ga\Q 2 ), for all three sets of DAs. Experimentally this ratio, i.e., G a{Q 2 ) I 'Ga (Q 2 ) 
should be close to 1. From these figures we see that the prediction for this ratio by LCSR 
is quite close to 1 at f3 = —1.4 for the first (third) set of DAs and for the choice of the 
mass m A = 1.0 GeV (m A = 1.2 GeV). The ratio R(Q 2 ) = G A (Q 2 )/G d A (Q 2 ) is close to 1 
for the second set of DAs at (3 = —1 for the choice of the mass itla = 1.2 GeV . But, as 
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has already been noted, (3 should be different from —1. Therefore we can conclude that the 
dipole form of the Ga{Q 2 ) describes well only the first set of DAs at (3 = —1. Our results 
on Ga{Q 2 ) are in close agreement with the chiral quark model [24] predictions and lattice 
results [25J at (3 = —1.4. The results for Ga{Q 2 ) for other values of (3 depart considerably 
from the lattice results. 

It should be remembered that all available experimental data exists only at low Q 2 
region. But unfortunately, our approach can give reliable prediction about the form factors 
only at high Q 2 > 2 GeV 2 region. Therefore direct comparison of our theoretical prediction 
on experimental result is impossible. Recently proposed experiment Minerva [26] would 
provide a precise determination of Ga at Q 2 < 2 GeV 2 and Q 2 > 2 GeV 2 . Another 
experiment using e + p — > v + n has been planned at JLAB [27], covering the range 
Q 2 = 1 4- 3 GeV 2 . When the result of these experiments will become available, it will be 
possible to compare the theoretical and experimental results. As the final remark we would 
like to note that, for a more reliable prediction of the form factors, it is necessary to take 
into account the radiative a s corrections and the distribution amplitudes with four, five 
particles. 

In conclusion, in this work we calculate the form factors of nucleons induced by isovector 
and isoscalar axial-vector currents in LCSR method using the general form of the nucleon 
interpolating current. In our calculations we use three different sets. The dependence of the 
form factor Ga{Q 2 ) on Q 2 is compared with its dipole form parametrization which follows 
from the analysis of neutrino experiments. Our analysis show the results for Ga{Q 2 ) in 
LCSR, for the first set of DAs at (3 — —1.4, well describes the dipole form of Ga{Q 2 )- 
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Figure captions 

Fig. (1) The dependence of the nucleon form factor Ga{Q 2 ) on Q 2 at M 2 = 2 GeV 2 and 
s = 2.25 GeV 2 , at four different values of (3: (5 — -5, (5 = -1.4 (5 — -1 and /? = 5, for 
the first set of DAs. 

Fig. (2) The same as in Fig. (1), but for the form factor Gp{Q 2 ). 
Fig. (3) The same as in Fig. (1), but for the form factor G S A (Q 2 ). 
Fig. (4) The same as in Fig. (1), but for the form factor G S P (Q 2 ). 

Fig. (5) The dependence of the ratio R = G a{Q 2 ) / 'G d A (Q 2 ) on Q 2 at three fixed val- 
ues of (3, for the first set of DAs at = 1.0 GeV. 

Fig. (6) The same as in Fig. (5), but for m A = 1.2 GeV. 

Fig. (7) The same as in Fig. (5), but for the second set of DAs. 

Fig. (8) The same as in Fig. (6), but for the second set of DAs. 

Fig. (9) The same as in Fig. (5), but for the third set of DAs. 

Fig. (10) The same as in Fig. (6), but for the third set of DAs. 
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